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We present direct numerical simulations (DNS) of the mixing of the passive scalar at modest 
Reynolds numbers (10 < R\ < 42) and Schmidt numbers larger than unity (2 < Sc < 32). The 
simulations resolve below the Batchelor scale up to a factor of four. The advecting turbulence is 
homogeneous and isotropic, and maintained stationary by stochastic forcing at low wavenumbers. 
The passive scalar is rendered stationary by a mean scalar gradient in one direction. The relation 
, between geometrical and statistical properties of scalar field and its gradients is examined. The 

Reynolds numbers and Schmidt numbers are not large enough for either the Kolmogorov scaling or 
the Batchelor scaling to develop and, not surprisingly, we find no fractal scaling of scalar level sets, 
£Nj ' or isosurfaces, in the intermediate viscous range. The area-to-volume ratio of isosurfaces reflects 

the nearly Gaussian statistics of the scalar fluctuations. The scalar flux across the isosurfaces, 
which is determined by the conditional probability density function (PDF) of the scalar gradient 
magnitude, has a stretched exponential distribution towards the tails. The PDF of the scalar 
dissipation departs distinctly, for both small and large amplitudes, from the lognormal distribution 
for all cases considered. The joint statistics of the scalar and its dissipation rate, and the mean 
conditional moment of the scalar dissipation, are studied as well. We examine the effects of coarse- 
graining on the probability density to simulate the effects of poor probe-resolution in measurements. 



(N 



Q 

U 
d 



o 
in 



X 



PACS numbers: 47.27.Gs, 47.53. +n, 02.70.Hm 



I. INTRODUCTION 



The mixing of passive scalars in the presence of turbulent motion is a subject of great theoretical and practical 
I . interest. Examples of its ubiquitous applications arise in reacting flows and combustion, mixing of salt and plankton 
' in oceans and of pollutants in the atmosphere, as well as mixing of biological substances. Much experimental effort 
£N| . has thus been invested on mapping scalar fields in three dimensions, usually by reconstructing them from several 
planar cuts taken in quick succession in time; see Refs. 1 to 6 for work on turbulent liquid or air jets and Refs. 7 and 
8 for work on diffusion flames. The Schmidt number Sc — v/n where v is the kinematic viscosity of the advecting 
fluid and n is the scalar diffusivity, plays an important role in determining the nature of mixing [9|. For the following, 
our focus will be on Schmidt numbers larger than unity. 

The finest scales of the scalar are generally thought to be of the order of the Batchelor scale defined as t\b = T)/ ^/~Sc 
where rj is the Kolmogorov scale given by r/ — i/ 3 / 4 / (e); e is the local energy dissipation rate of turbulence, and the 
angular brackets indicate a suitable average. For a given flow, the Batchelor scale becomes smaller for larger Sc. 
Thus, the required resolution for large 5'c, common in liquid flows, is hard to attain experimentally |Toj . 

With modern computers, it is possible to resolve scalar mixing in turbulence at high Schmidt numbers [ill IT^ . IT3I 
0,^3- Schumacher et al. ll| showed that it is sometimes necessary to resolve scales that are finer than r\B- Such fine 
resolutions are desired, for instance, for obtaining reliable results on the extreme fluctuations (both small and large) 
of scalar gradients and scalar dissipation rate. We resolve scales finer than the Batchelor scale by a factor of four. 
This fine resolution limits the ranges of Sc and R\ examined, where R\ is the Taylor microscale Reynolds number. 
Nevertheless, we can study trends with Sc and R\ in limited ranges of these parameters, with full three-dimensional 
field resolved extremely well. In particular, we consider relations between geometric and statistical properties of 
passive scalars and their spatial variations with respect to Sc for a given turbulent flow, as well as those with the 
Reynolds number for a given Sc. 

The paper is organized as follows. Section II briefly reports the numerical details and parameters. We consider the 
box-counting properties and area-to-volume ratio of the passive scalar isosurfaces in Sec. Ill, while Sec. IV considers 
the properties of the instantaneous scalar flux across isosurfaces. Section V presents an analysis of the statistics of 
scalar dissipation including conditional statistics. Conclusions and summary are presented in the final section. 



II. NUMERICAL SIMULATIONS 



We solve simultaneously the Navier-Stokes equations for an incompressible flow u(x, t) and the advection-diffusion 
equation for the passive scalar field #(x, t) in box-type three-dimensional turbulence with the classical pseudo-spectral 
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0.307 0.345 


0.372 


0.306 
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23.68 


11.84 5.92 


23.68 11.84 


5.92 


2.82 


Le 


0.707 


0.587 0.477 0.562 0.447 0.357 0.289 



TABLE I: List of parameters for the present DNS runs. N is the number of grid points on each side of the computational cube, 
v is the kinematic viscosity, (e) is the mean energy dissipation rate, u rms = (u 2 ) 1 ^ 2 is the root-mean-square velocity fluctuation 
in x direction, k m = \/2N/3 is the maximum wavenumber allowed in the simulations, R\ — yjl§/(v{t)) (u 2 ) is the Taylor 
microscale Reynolds number, L is the integral scale of turbulence. Its magnitude is to be compared with the linear dimension 
of the box, which is 2n in all cases. T av is the statistical averaging time, to be compared here with the large-eddy turnover 
time defined as Te = (u 2 )/2(e); Sc is the Schmidt number, 6> rms = {O 2 ) 1 ^ 2 is the root-mean-square passive scalar fluctuation, 
(te) is the mean scalar dissipation rate and Le is the integral scale of the passive scalar. Both integral scales are calculated as 
in Ref. 44. 



method that uses fast Fourier transformations and a 2/3 de-aliasing. The boundary conditions for the velocity and 
scalar field are both periodic. The simulation domain has a length of 2ir in each coordinate direction and is resolved 
by an equidistant grid with N 3 points. The equations are 

^ + (u-V)u = -Vp + i/V 2 u + f, (1) 
Vu = 0, (2) 
— + (u-V)0 = nV 2 6-u y G. (3) 

Here, p(x , t) is the kinematic jpressure field. The turbulence is sustained by a random force density f (x, t) (see Vedula 
et al. [llj and Schumacher yjj for details) and the scalar gradient is constant in the y direction, i.e. G = (0, G, 0) 
with G = 1; u y is the turbulent velocity component in the direction y. The equations are integrated in time by 
a second-order predictor-corrector scheme. Parameters of the simulations and some statistical results are listed in 
Table 1. The spectral resolution exceeds the usual criterion of k max r\B > 1.5 by a factor of about 4. A detailed 
comparison between these high-resolution data and those with the nominal resolution employed in most DNS studies 
was presented in Schumacher et al. |l6j | . Their conclusion was that super-fine resolution would be needed to capture 
specific features of the scalar dissipation rate 

e e (x,t) = k(W(x,£)) 2 , (4) 

such as the sheet-like structures in which the extreme events are found. We shall exploit the present fine resolution 
for just such purposes. 



III. SCALAR ISOSURFACES 

A. Box-counting dimension 

It is well-known that, for Schmidt numbers exceeding unity, the passive scalar is filamented in structure and 
mixed by smooth flow in a continual stretch-twist-fold scenario |18| . An open point is whether the scalar contours 
accompanying such dynamics, when diffusion and external driving are included, generate fractal sets. Neither the 
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FIG. 1: Two-dimensional instantaneous slice through the fluctuating passive scalar field #(x,£) for the level of O.S8 rms (i.e., 
X — 0.8). The side-length of the box (Run 2) is about 90r] and resolved with N = 1024 in each coordinate direction. R\ = 24, 
Sc = 32. 
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FIG. 2: Three-dimensional box-counting data for scalar level sets Lx as defined in Q for Sc = 32 and for R\ — 10 (left 
column) and R\ — 24 (right column). The values of X are: asterisks, 0.8; circles, 1.6; squares, 2.4; triangles, 3.2. The vertical 
dashed lines in the upper panels show the corresponding values of t]b- Short solid lines mark the algebraic slope with 5~ 2 in 
the log- log plot. All the isosurfaces have a region with a common exponent of 2 at the smallest scales. The lower panels show 
the local slope d log(Ng)/d log(S) as a function of S. Symbols are as in the upper figures. 
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inertial-convective range nor the viscous-convective Batchelor range obtains for the present moderate values of R\ 
and Sc, so it would be surprising if the fractal scaling is found. Nevertheless, the issue is worth a brief examination 
as a prelude to the study of area-to- volume ratio of scalar interfaces that follows in Sec. Ill B. If the surface is indeed 
a fractal, it is easy to write down the surface area in terms of its dimension and the cut-off scales terminating the 
scaling behavior |19|. 

Let us define the scalar level sets as 

L x = {x : X - SX < 0(x, t)/0 rms < X + SX} , (5) 

where X is the ratio of the chosen amplitude threshold to the root-mean-square value of the scalar, 9 rms . Operationally, 
5X is taken to be 2.5% of rms for each (R\,Sc). Figure 1 shows a level set of the scalar field around a value of 
0.88 rms for Sc — 32 and R\ = 24. Similar pictures with different thresholds show that levels of high intensity are 
embedded within iso-levels of lower intensity. The iso- levels become increasingly disconnected with increasing intensity 
and become isolated islands embedded within an iso-contour of lower magnitude. This observation is consistent with 
mixing studies in a jet by Villermaux and Innocenti In the box-counting procedure described, for example, in 
Falconer [2(J], one counts the number of cubes, Ng, of side-length 6 that cover the level set Lx completely. This 
procedure is repeated for various values of S. The box-counting dimension D of a (fractal) level set is defined then as 

D = -lim^. (6) 

<5^0 logo 

In practice, it is clear that one cannot apply the limit S — ► and one seeks an algebraic scaling, N$ ~ 5~ D , in a finite 
range. 

The upper two panels of Fig. 2 show the log-log plots resulting from box-counting for data sets of Runs 1 and 2 at 
Sc — 32. A constant slope with exponent D would give the box-counting dimension of the scalar level set. However, 
the local slope given by 

dfog^ 
dlogo 

varies continuously, as shown in the lower panels of Fig. 2. Though our simulations cover three orders of magnitude of 
scales, it is clear that no power-law scaling exists for reasons mentioned at the beginning of this section. The smallest 
few values of 5 give a local slope of 2, for all X, which means that the level sets at those scales become smooth sheets. 
We interpret (the slightly) smaller slopes for intermediate scales of 5 as arising from the disconnected shape of the 
iso-level sets. Indeed, this patchiness increases with the level X thus giving a smaller magnitude of the local slope. 

The conclusion from this study is that fractal behavior of the isolevel sets is not obtained in the range of R\ and Sc 
considered here. We want to emphasize, however, that this finding has no bearing on the past experimental work as 
reported in Refs. 3, 4, 19 and 21. The differences such as may exist in these experimental studies require a separate 
discussion which lies outside the scope of this paper. We expect that the fractal property in the inertial-convective 
range requires high Reynolds numbers while that in the viscous-convective range requires high values of Schmidt and 
Reynolds numbers. 



B. Area-to-volume ratio 



Since we do not have fractal scaling for the ranges of parameters considered here, it is necessary to compute the 
area of isosurfaces directly. In reacting flows, a quantity of some importance is the ratio of the isosurface area content 
to the volume which encloses it 22] . This ratio is numerically calculated as follows. We first define the set 

L x ={x:0(x,i)/0 rms >X}. (8) 

The set is computed for each snapshot separately and the results are averaged with respect to time subsequently. We 
mark all grid points that belong to Lx by an indicator function 1^. = 1. For all other points, ly^ = 0. Here i, j, 
and k denote the indices of the grid vertices running from 1 to N. All boundary points of the set Lx have a number 
of nearest neighbors, Nij k , that ranges from 1 to 5; the inner points of Lx have 6 nearest neighbors. The number of 
inner points is denoted by Nijk- The area-to- volume ratio Tx of the set Lx is given by 
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FIG. 3: Left column: mean of the area-to-volume ratio (IV) as a function of the isolevel X in a log-log plot. Upper panel 
shows data for Run 2 at three different Schmidt numbers. Lower panel shows the quantity at Sc = 32 for Runs 2 and 3. Right 
column: the corresponding mean thickness scale (£x) as derived from Eq. (lli)l . 




FIG. 4: Sketch of the scalar gradient flux calculation across an passive scalar isosurface X = Xo. The illustration is for a 
two-dimensional case (but the data presented in the next figure are for the surface in three dimensions). Black points lie on the 
boundary whereas the gray points are the interior of Lx- The inner points are the volume that is used for Tx (see Eq. @). 



with grid spacing Ax, with Nij^ ^ 0. The left column of Fig. 3 shows the average of the area-to-volume ratio, with 
(Tx) taken over several snapshots. Three trends are clear: first, the ratio grows with increasing X, which indicates 
a higher degree of disconnection of the isosurfaces. Second, the ratio is larger for larger Schmidt numbers. This is 
because mixing then occurs over ever finer scales whose contours are stretched and folded by the flow. Finally, a 
larger Reynolds number enables the scalar to be stirred more efficiently on all scales, yielding larger area-to-volumc 
ratios with increasing Reynolds number. 

We relate this geometric measure to statistical properties of the passive scalar fluctuations following ideas that were 
used in the slightly different context of interfacial dynamics for a turbulent jet [23. As is shown later (see Sec. V 
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FIG. 5: The probability density function of the magnitude of scalar gradient, |VX|, taken at iso-contour levels X — Xo- Upper 
left: The data are for Run 3, Sc — 32, for three different levels: X — 0.5 (open circle), X — 1.0 (asterisk), X — 1.5 (triangle). 
Upper right: The data are for Run 2, Sc = 32, for three different levels: X = 0.5 (open circle), X = 1.0 (asterisk), and X = 1.5 
(triangle). Lower left: Schmidt number dependence of the flux. Data are for Run 2 at X = 1: Sc — 2 (circle), Sc = 8 (asterisk), 
and Sc — 32 (square). Lower right: Reynolds number dependence of the flux for Sc = 32 and for X — 1, with R\ = 24, Run 2 
(square) and R\ = 42, Run 3 (circle). The vertical arrows (dashed for Run 2 and solid for Run 3) indicate X/tjb with X = 1. 



B), the scalar fluctuations are close to the Gaussian probability density function, p(X) = \/\/2tt exp(— X 2 /2), with a 
slight sub-Gaussian behavior in the far tails. The average of the ratio (Fx) should therefore be related to p(X). For 
X = Xq we set 

, r ue > exp(-X 2 /2)(LYo exp(-X 2 /2)dX 

(l Xo) { l Xo) _ ^ exp( _ y2/2) dy erW x/2)] ' 1 Uj 

where erf(x) = (2/y / 7r) J* exp(— y 2 )dy is the error function. For the following, we discuss results for values X > 0, 
but everything holds for Xq < just the same way. The additional scale (ix ) on the left hand side of Eq. (10) has 
to appear for dimensional reasons since the right hand side is dimensionless. 

The scale (£x ) is related to the surface-to- volume ratio T, and is small when the surface-to- volume ratio is large. 
Since T becomes larger for "rougher" surfaces, the scale is a characteristic measure of the transition between smooth 
and rough isosurfaces, perhaps describing the average connectedness and regularity of the isosurface. We may thus 
call it the transition scale. It is a global average measure as well as T. Its trends with Sc and R\ are shown in 
the right column of Fig. 3. The straight horizontal parts of (£x) for X < 2 at fixed Sc and R\ are the "geometric 
fingerprints" of the Gaussianity of the scalar fluctuations. Finally, we want to note that the relation l|10fl holds only 
for isosurfaces that are not fractal. 



7 




20 40 60 80 100 120 



z =V< e e> 

FIG. 6: Log-linear plot of the probability density function of the scalar dissipation rate, normalized to the mean value. Data 
are for Run 1 and Run 2 at Sc = 32. Fits to the data for z > 10 with the exponential term of l)13|l are also plotted and the 
corresponding exponents a are shown. The dashed line is the optimum of a least square fit resulting in a = 0.86. 

IV. FLUX ACROSS SCALAR ISOSURFACES 

One reason for the interest in the area of an isosurface is the flux across it. For example, in non-premixed turbulent 
combustion the isosurface of the stochiometric mixture fraction plays an important role and the variation of gradients 
across it are important for closure models [lfl I22I l24j . This can be calculated by integrating the flux across all 
infinitesimal elements of the surface. Since the flux is proportional to the product of the area and the gradient of 
the concentration across it, it is necessary to know the gradient across each infinitesimal element of the surface. The 
gradient can itself be a highly fluctuating variable. 

The present data allow the computation of gradients and fluxes across infinitesimal elements of isosurfaces. The 
(differential) flux across an infinitesimal area element dA of an isosurface is given by 

d$(X ) = -VX\ x=Xo -ndA, (11) 

where n is the outward unit normal (see Fig. 4 for Xq > 0). 

The probability density function (PDF) of d$(Xo) is then fully determined by the PDF of scalar gradient magnitude 
at the isolevel X — Xq 

p(d<S>(X )) = p(-\X7X\ x=Xo dA) . (12) 

Thus, the quantity of interest is a conditional PDF p(|VX| \X — X ) which is calculated as follows. First, as in the 
preceding section, we mark all the grid points for which 0(x, t) > X9 rrns as active points; their grid sites are marked 
with Iij k = 1; for all other points, Iij,k = 0. Gradients are calculated on dL x afterwards. Sparser grid resolutions 
would require a trilinear interpolation between the grid cells |25| . 

The result is plotted in Fig. 5. A stretched exponential behavior is observed for a large range of the PDF, with 
the tails decaying more rapidly for larger X. The PDFs suggests also that there is a considerable local variation of 
the flux across interfaces X = Xq. The fluctuations increase significantly with increase in Sc as well as R\. For both 
Reynolds numbers, one can see low-probability instances whose gradients are of the order of 9 rms /rjB- This particular 
gradient value is indicated in the lower right figure by the two vertical arrows. 



8 




-6 -4 -2 2 4 
z=[ln(e 9 )-<ln(e e )>]/o 




z=[ln(8 e )-<ln(e e )>]/a 

FIG. 7: Upper panel: Probability density function of z — (log(eg) — (log(ee)))/<7 with a = -J ((log(e9)) 2 } — (log(eg)) 2 . Data 
for Run 2, with Sc = 32, are compared with the corresponding lognormal distribution having the same mean and standard 
deviation. Data in each panel are plotted as a solid line and the lognormal distribution as a dotted line. The inset shows the 
same data in a log-linear plot that highlights the tails. Deviations from lognormality are quite visible. Lower panel: Log-linear 
plots of the same type of data for three different Reynolds numbers and Sc = 32. 

V. STATISTICS OF SCALAR DISSIPATION RATE 

A. The PDF of the scalar dissipation rate 

The scalar dissipation rate plays a central role in determining turbulent mixing. In combustion, for example, its 
properties enter fast chemistry models [2(j , modelling of thin reactive layers called flamelets (2?| which are embedded 
and advected in the turbulent flow, and conditional-moment closures for the evolution of the mixture fraction in 
non-premixed cases [28|. 

With p(|V0|) ~ exp(— (C| V6\ a ) for large |V#|, as was seen in the last section, we expect that the tails of the PDF 
of eg will have stretched exponential behavior given by 

p(e e > (e e )) = -%exp(-C 2 e° /2 ) , (13) 

where C\ and C2 are constants. Such statistics were derived analytically for scalar advection in smooth and white-in- 
time flows for the limit Pe — > 00 and a = 2/3 was found |2^,|33|. In Fig. 6 we show the PDF of the scalar dissipation 
rate in a log-lin plot for Runs 1 and 2; Sc=32. In addition, we fit the exponential part of the formula (|13l) to the data 
for the range of eg > 10(ee). The cases a = 1 (pure exponential) and a — 2/3 (corresponding to the theory in Refs. 
29 and 30) bound both sets of data from below and above, respectively, with the optimal least-square fit (dashed line) 
corresponding to a = 0.86. 



9 




z=[ln(E e Mln(E e )>]/o 



FIG. 8: Log-linear plot of the PDF of the scalar dissipation rate at Sc — 32 for Run 2. As in Fig. 7, solid line are the original 
data at a resolution oi N = 1024 and the dotted line is for the corresponding lognormal distribution. The analysis was repeated 
for the coarse-grained scalar dissipation field, and the dashed line shows the result of coarse-graining over 8 3 cubes, giving an 
effective resolution of N = 128; the dashed-dotted curve is for 4 3 cubes with an effective resolution of N = 256. 

This form of the PDF departs clearly from log-normality. Figure 7 (upper panel) shows the PDF of eg for Run 2 
(see Table 1). The linear scale reveals that the lognormal approximation is not particularly good near the core. The 
logarithmic scale (inset) shows that the tails depart from lognormality as well: low dissipation values show a fatter 
tail and high dissipation values fall short of the lognormal. These departures from lognormality are robust for all 
Reynolds and Schmidt numbers considered here: this is demonstrated in the lower panel of Fig. 7 where we plot the 
data for Sc = 32 at the three Reynolds numbers i?A = 10, 24, and 42. 

There has been much discussion in the literature about the lognormality or otherwise of scalar dissipation, and 
we might cite an early paper Sreenivasan and Antonia [3l| and a recent one by Su and Clemens to illustrate 
experimental results and the discussions accompanying them. Often, measured deviations from lognormality are 
attributed to poor measurement resolution or to low Reynolds numbers or to noise in the experimental system. 

In order to shed some light on how the PDFs depend on probe resolution in an experiment, we plot in Fig. 8 
the PDFs that result when the statistics are taken over successively coarse-grained scalar dissipation fields. The 
coarse-graining scale does not produce a qualitatively different result on the right side of the distribution (z > 0), 
for which the PDF remains sub-Gaussian for all cases. For negative z, there is probably a particular coarse-graining 
that may collapse the data on to the lognormal curve, but this coincidence would be accidental at best. On the effect 
of Reynolds number, we can only comment that the behavior is similar for both Reynolds numbers considered here. 
The double-precision calculations throughout this work keep the noise effects under control. We thus conclude that, 
while the PDF is roughly lognormal, the deviations from it are considerable. 

B. The joint statistics of scalar and scalar dissipation 

In Fig. 9 we show the joint PDF p(9, eg) where p(9, eg) d9 deg is the probability of having the scalar dissipation rate 
between eg and eg + deg and the scalar fluctuation between 9 and 9 + d9. The figure shows the joint PDF for two 
Reynolds numbers and Sc = 32. The upper row of Fig. 9 shows iso-contours of p(9, eg) which are equidistant in units 
of the logarithm to base 10. The joint PDF is wider with respect to the 6>-axis for the case of larger Reynolds number. 
This is due to the fact that the scalar field and the dissipation field acquire larger amplitudes with increasing R\ . The 
lower panels show the scalar PDF calculated by integrating p(9,eg) over eg. We observe an almost Gaussian shape 
of p(9) for both Reynolds numbers (as anticipated in Sec. III.B). Additionally, we have compared our results with 
the limiting shape for a scalar field PDF that was analytically predicted for advection in a random flow |32|. The 
comparison is good only for small amplitudes of 9. 

A nearly Gaussian PDF for the scalar fluctuations is in agreement with a number of other studies at Sc < 1, 
such as heated grid turbulence |3^| , isotropic turbulence with a mean scalar gradient 0> > a homogeneous shear 
flow with heating by a constant mean temperature gradient |36| . In our simulations, the fluctuations extend to 
approximately ±4.56*,. ms and the deviations in the tails are slightly sub-Gaussian. Exponential tails of the PDF, as 
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found by Jayesh and Warhaft [33, were not observed for our range of parameters. It should be mentioned that they 
were also predicted analytically for a white-in-time smooth Gaussian flow by Shraiman and Siggia 38]. Furthermore, 
exponential tails have been suggested to result from an aggregation process of elementary local straining motion 
events, and close-to-exponential tails have been observed in jets sufficiently far away from the nozzle exit [39(. 

Overholt and Pope yj] argued that the nearly Gaussian or slightly sub-Gaussian statistics may result from the 
limited size of the simulation box in relation to the integral scale. The ratio of our box size to the integral scale of 
the flow, 2ir/L, is about 6 to 8 (see Table 1); in terms of the integral scale of the passive scalar at the largest Schmidt 
number, the ratio 2ir/Lg varies from 13 to 21. This is significantly larger than the corresponding scale ratios in the 
experiment of Jayesh and Warhaft. The available evidence thus suggests that the relative size of the box size to the 
integral scale is not responsible for the shape of the PDF. A plausible argument can be made that the ratio L/Lg 
might be a more appropriate quantity. We have attempted to vary this ratio in the simulations, but the results are 
still not conclusive. 

The average of the scalar dissipation rate conditioned on the relation 0(x,t) = ip is denoted by (eg\9 — ip)- The 
conditional average is defined by means of the conditional probability density function p(eg\9), as 

(ee\6 = Tp) = / de e e eP (e e \e = V) • (14) 
Jo 

The conditional probability density function can be expressed via the joint probability density function of having the 
scalar dissipation between eg and eg + deg and the scalar between 9 and 9 + d0, so that 

p{ee\6 = ^)= p{0 = ^ ■ (15) 

Figure 10 shows the conditional moments for Sc — 2 and 32. It is well known that the conditional moment is 
constant |40j | if the scalar distribution is strictly Gaussian. As can be seen in Fig. 10 our numerical results confirm 
this prediction for scalar fluctuation levels of less than about three standard deviations. For larger scalar fluctuation 
levels, there is a considerable drop in the conditional expectation. This suggests a precipitous drop in eg for large 
values of 9, in contrast, say, to the behavior of the energy dissipation field e with respect to the velocity fluctuations 

m _ 

Experiments in inhomogeneous flows |42j show that the constancy of the conditional moment holds roughly on the 
centerline of turbulent jets and wakes. The behavior off-center is quite complex, and the boundary layer flows are 
intrinsically different. 



VI. SUMMARY AND CONCLUSIONS 



We have studied numerically the dynamics of passive scalars advected in homogeneous isotropic turbulence for three 
Schmidt numbers Sc > 1 and three Reynolds numbers, R\. These Reynolds and Schmidt numbers are moderate but 
comparable to those of many laboratory experiments on reacting flows. Usually, in such flows, neither the inertial 
range nor the viscous-convective range is well developed. The most important range of scales then belong to the so- 
called intermediate viscous range, which essentially pertain to the crossover between the inertial and viscous ranges 
of scales |l6l |43j. In this range, rough velocity filaments reach down into the viscous range below the Kolmogorov 
scale T) and stir the scalar. As reported in Ref. 16, the intermediate viscous range gets more extended with increasing 
Reynolds number. We have paid particular attention to the grid resolution of the simulations, which exceeds the 
nominal value by a factor of four, to ensure that the high intensity events of the scalars and their gradient fields are 
resolved well. 

The first quantity considered is the area-to-volume ratio of isoscalar surfaces because of its fundamental relevance 
to reacting flows. For the Reynolds and Schmidt numbers considered here, the fractal scaling of isoscalar level sets is 
not obtained and so cannot be used to estimate the area-to- volume ratio. This ratio, computed here directly, increases 
with an increase in the passive scalar isolevel, the Schmidt number, and the Reynolds number, respectively. In detail, 
these increasing trends depend on the statistics of the scalar itself, which is found here to be nearly Gaussian. 

We have defined a so-called transition scale (£x) that becomes smaller with increasing surface-to- volume ratio of the 
isosurface corresponding to the scalar magnitude X (normalized by the root-mean-square value). Quite understand- 
ably, the isosurface becomes "rougher" and more disconnected with increasing Sc and increasing Reynolds number, 
but is not dependent on X itself — except when X is large. The flat behavior of the scale (£x) for X < 2 can be 
considered a geometric confirmation of the nearly Gaussian probability density function p(X) of the scalar. For large 
X, there is a tendency for (£x) to decrease with X, this being attributable to the sub-Gaussian behavior of p(X) 
in the tail region. The flux across the interface is a strongly fluctuating quantity from one differential area element 
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FIG. 9: Upper row: Joint probability density function p(9/9 rms ,eg /{eg}) for the runs at R\ — 24 (left) and R\ = 42 (right), 
both for Sc — 32. The iso-contour levels are equidistant in units of the logarithm to base 10, varying from —6 to 0.5 in steps 
of 0.5. The joint densities for large negative values of log(ee) (i.e., scalar dissipation rates smaller than 10~ 3 (eg}) are truncated 
for reasons of poor convergence. Lower row: scalar PDF for R\ = 24 and Sc values indicated in the legend (left). The thin 
dashed line is a theoretical prediction for the scalar PDF, p(X) = C/(l + kX 2 ) 1+1/l2k , for advection in a random flow with fit 
parameters C = 0.4, k = 0.08, X being 9/9 rma [32^1 . The thin dotted line is a Gaussian distribution. Scalar PDF at Sc = 32 
for two different Reynolds numbers as indicated in the legend are shown in the right panel. Both thin lines are as in the lower 
left panel. 

to another, and its distribution is closely exponential. The fluctuations increase rapidly with increase in Sc as well 
as R\. Though with low probability, gradients of the order of magnitude rms /rjB do occur in the scalar field. This 
clearly suggests the existence of sharp fronts accompanied by large jumps in the scalar concentration. 

The scalar dissipation rate departs from lognormality more than previous measurements suggest. The reasons for 
close correspondence in previous data are not entirely clear, but it must be pointed out that very few data sets 
exist in which all three scalar gradient components that enter the scalar dissipation rate definition (see (Eq. 0J) 
were measured with adequate resolution and without the use of Taylor's hypothesis. To consider if the closeness to 
lognormality observed previously was caused by poor resolution in measurements, we coarse-grained the data and 
replotted them. The coarse-grained data also display significant departures from lognormality. We thus conclude that 
the reasons for observed departures are more complex; they may also emphasize the lack of strict universality in the 
dissipation statistics. 

The high spectral resolution in the present DNS allowed the evaluation of a conditional moment analysis as well. 
We have been able to confirm theoretical predictions that follow from nearly Gaussian fluctuations of the passive 
scalar field for 9 < 39 rms , say. 

In summary, we have quantified various properties of the scalar isosurfaces for conditions that are typical of many 
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FIG. 10: Mean conditional moment of the scalar dissipation rate for Schmidt numbers Sc = 2 and 32. Data are for Run 2 in 
Table 1. 



laboratory experiments in reacting flows. Beside the connections between statistical and geometrical properties of 
scalars studied here, it is interesting to shed more light on the formation of characteristic scalar structures with respect 
to time, and to relate them to the underlying local flow properties for different Reynolds numbers. This can be a 
useful way of incorporating Reynolds number effects on mixing in the Batchelor regime and is part of the future work. 
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